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: ABSTRACT 

^ ■ We have mapped the high-redshift (z = 0.546) cluster CI 0016+16 with the Ryle 

' Telescope at 15 GHz. The Sunyaev-Zel'dovich decrement is clearly detected, and 

\ resolved. We combine our data with an X-ray image from ROSAT, and a gas temper- 



ature from ASCA to estimate the Hubble Constant Ha ~ 69_[g km s Mpc for an 



^ ■ = 1.0 cosmology or //() = 84+ j9 km s"' Mpc"' for =0.3 and i^A = 0.7 

Q ■ Keywords: cosmicmicrowavebackground-galaxies:clusters:individual:C10016+16 
■ 

' - galaxies: intergalactic medium - radiative transfer - X-rays: sources 



1 INTRODUCTION 



The cluster CI 0016-1-16 lies at a redshift of 0.546 ( Pressler & Gunn 1992| ), and in X-rays is 



one of the most luminous clusters known ( JHenry et al. 1992| ). Despite its high redshift it has a 



predominance of red galaxies (i.e. it defies the Butcher-Oemler effect ( [Butcher & Oemler 1984| )). 



There is also evidence that it is undergoing a merger; the galaxy distribution is distinctly bimodal, 
and a map of the total mass derived from weak tensing of background galaxies ( [Smail et al. 1995| ) 



shows sever al peaks, although the X-ray emission appears somewhat smoother (B irkinshaw & 
Hughes 1994). ' 

CI 0016-1-16 was also one of the first clusters in which the Sunyaev-Zel'dovich (SZ, Sunyaev & 
Zel'dovich (|1972[)) effect was detected ([Birkinshaw et al. 1984]). It has been subsequently mapped 



with the OVRO/BIMA array at 28-30 GHz ( [Carlstrom et al. 1996D and detected at 2.1 mm and 



1.2 mm with the Diabolo experiment at the IRAM-30 m telescope ( pesert et al. 19981 ). Here we 
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2 Keith Grainge et al. 

present Ryle Telescope (RT) observations of the SZ effect in CI 0016+16 and provide an image of 
its two-dimensional structure. We then combine this data with the ROSAT PSPC image and ASCA 
temperature data to calculate the Hubble constant. 



2 SZ OBSERVATIONS WITH THE RYLE TELESCOPE 

The RT is an east-west aperture synthesis telescope specifically designed to map low-surface- 
brightness structures in the microwave background ( pones 199T1 ). Five of its eight antennas are 
used in close-packed configurations to provide the high temperature sensitivity needed to map 
the SZ effect. The system temperature at 15 GHz is 60 K, giving a flux sensitivity in 12 h of 
200 //Jy beam ^ The RT correlator provides 35 spectral channels of approximately 10 MHz band- 
width each (7 synthesised channels in each of 5 IF sub-bands). 



We have made many detections of the SZ effect with the RT (P^ones et al. 1993|; Grainge et 



I. _ • . I 



al. 1993; [Jones 1995| ; [Saunders 1995D ; we will first describe in some detail the general observing 



procedure, and then discuss the CI 0016-1-16 observations. 

2.1 General observing procedure 

Clusters are observed in runs of 12 hours with the telescope in one of the compact configurations 
shown in table [l|. In the case of fields of declination less than 48° with the RT in configuration Ca 
or Cb, the observation is curtailed when none of the baselines of the compact array remain unshad- 
owed. Due to the east-west nature of the RT, this results in a gap along the v-axis in the aperture 
plane coverage. To fill this region we also observe in either Cc array (for 24° < dec < 48°) or Cd 
array (for 16° < dec < 24°). The pointing centre chosen for each observation is the position of the 
peak of the X-ray emission from the cluster. For each observation, a bright (> 0.2 Jy) nearby point 
source is chosen as a phase calibrator and is observed for either five or ten 32-second samples 
(depending on its brightness), for every 30 samples on-source. A primary calibrator, either 3C48 
or 3C286, is observed either before or after a run. 

2.2 Editing the data 

The Postmortem package ( [Titterington 199 ip is used for editing and calibrating of RT data. 
The performance of the RF system is monitored by the use of a modulated noise source at each 
feed-horn. Corrections to the data and their weights are applied to take account of variations in 
system noise, mainly a result of weather conditions. Exceptionally bad data are discarded, along 
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The SZ effect in CI 001 6+16 3 

with that during winds in excess of 20 knots (marginal pointing) and during slewing and settling 
while changing pointing centres, and for any baseline actually, or close to being, shadowed. 



2.3 Primary calibration 

The absolute flux scale of the RT is fixed on 3C48 and 3C286. Flux densities of 1.70 and 3.50 Jy 



(Stoke's I — Q) are assumed for these respectively dBaars et al. 1977D . We calculate a set of com 



plex gains to apply to every sub-band and channel of each antenna from the primary calibrator 
observations. A least-squares algorithm is used to derive the gains for each antenna from the data 
for each baseline. This means that data from many baselines are used to calibrate each antenna, 
and increases the signal-to-noise of the calibration. The amplitude scale is observed not to drift 
significantly over the timescale of individual observations, and so is calibrated with a single cor- 
rection for each observation. This primary calibration is then applied to the data from the main 
observation. Monitoring by the VLA has detected low level variability in both of these sources 
on long time-scales (JVLA calibration manual 1996[), and we apply a separate correction for these 



changes. 

2.4 Phase calibration 

Antenna-based phase corrections are calculated for each of the 5- or 10-sample secondary calibra- 
tor observations, and interpolated to correct the source phases. This removes phase drifts due to 
the atmosphere and local oscillator distribution, which can be up to 360° during a 12-hour run for 
the longer baselines of the compact array. Since phase differences between channels in the same 
sub-band are relatively constant, and are calibrated out using the flux calibrator observations, we 
merge the visibilities of the 7 channels for each sub-band and smooth over all 5 or 10 samples 
within the particular calibrator pointing in order to increase the signal-to-noise ratio of the calibra- 
tion. We do not merge the secondary calibrator visibilities over the sub-bands since we have found 
that during a 12-hour run the different sub-bands can drift by up to 5°. Closure phase errors are 
calculated when the antenna solutions are found, and are almost always consistent with the thermal 
noise. The data are flagged where this is not the case. In order to check the secondary calibration 
the visibilities are smoothed, their mean amplitudes and phases are examined and the data for one 
representative channel are plotted out. We pick out and flag anomalous mean amplitudes, phases 
and their rms, and look for amplitude or phase jumps caused by interference or possible hardware 
faults. 
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2.5 Smoothing, weighting and interference clips 

The data are smoothed to an integration time appropriate to the resolution of the array: 320 s for 
the most compact arrays Ca, Cb and Cc, and 96 s for array Cd. These smoothing periods were 
chosen to ensure that Ra, the fractional decrease in the peak response to a point source is close to 
unity. Ra is given by (see e.g. Thompson ( |1986[ )) 

K^i-li'-^^X. (1) 



3 V 

where (Og is the rotation rate of the Earth, Xq is the averaging time, Qb is the FWHM of the synthe- 
sised beam, and r is the distance of the point source from the observation centre. For the smoothing 
periods used, a point source 6 arcminutes from the pointing centre (the position of the first null 
of the RT primary beam response) will have Ra > 0.99 on the longest baseline available in the 
respective configurations. Any low-level interference is removed by a simple clip at 0.25 Jy before 
smoothing, and at 0.079 Jy after smoothing over 10 samples, or at 0.144 Jy after smoothing over 
3 samples. This corresponds to approximately 3 standard deviations. The visibilities are weighted 
by their inverse variances, calculated from the product of the system temperatures of the two an- 
tennas and are then written out in FITS format to be transferred to AlPS. As a final check of our 
calibration procedure, each day's data are mapped individually using the entire av-range and natu- 
ral weighting, in order to give maximum signal to noise. We examine these maps for artifacts due 
to poor calibration and interference which have not been removed. 

2.6 Data analysis 

Before a map of the SZ effect can be made, the effects of discrete sources on the data have to be 
removed. The clusters we have studied have been pre-selected to be radio-quiet, but may contain 
several radio sources with flux densities of up to a few mJy. Rather than attempt to identify and 
determine the fluxes of these sources in the map plane through CLEANing the long baseline data 
as has been our practice previously dGrainge et al. 1993| ; [Tones et al. 1993[ ), we now fit directly to 



the visibilities in the aperture plane using the program Fluxfitter. This procedure is explained 
in greater depth in Grainger et al. ( [200 ID . We subtract the sources from the original Mv-data using 
point-source models, implemented by the AlPS task UVSUB. The procedure that we adopt for 
removing sources which are not point sources on our longest baselines, but are not extended on 
the scale of the SZ decrement, is detailed in Grainge et al ( |1996| ). If we find a bright source which is 



extended on the scale of the SZ effect, we cannot be confident of subtracting the source accurately 
and observations of the field have to be abandoned. 
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The SZ effect in CI 001 6+16 5 

The SZ decrement is typically several arc-minutes in extent and so is resolved out by the long- 
baseline antenna pairs, and to a good approximation is only detected by the short baselines. We 
therefore either map only the short baseline source-subtracted data (typically 0-1 kX), or apply a 
Mv-taper to look for an SZ decrement close to the X-ray cluster centre. The quality of the resulting 
image is determined by both the structure of the SZ decrement and the aperture-plane coverage of 



the observ ation. For high-declination objects, e.g. Abell 2218 (JJones 1995[) or Abell 773 (G rainge 



I I 

et al. 1993), the aperture-plane coverage is basically a circular ring at the radius of the shortest 

baseline, with very little SZ signal at the second-shortest baseline. An image made from this will 
show the position of the centre of the decrement and reveal any significant departures from circular 
symmetry, but does not strongly constrain the shape of the SZ profile. Additional information, 
such as X-ray imaging and assumptions about the physics of the gas, must be used to fully model 
the gas structure. In a low-declination object such as CI 0016-1-16 however, the projection of the 
baselines results in a much better-filled aperture plane. Using two different array configurations 
optimised for the north-south and east-west parts of the aperture plane results in an image in which 
significant structure can be seen. Since the SZ signals generally have only moderate signal-to-noise 
(~ 5a), CLEANing to remove the point-spread-function has to be performed with great care. With 
the poor aperture-plane coverage used for the mapping, the inner part of the synthesised beam is 
poorly represented by the Gaussian restoring beam, with the result that the map's noise level will 
be artificially reduced, and that significant artifacts can occur. The best results are obtained by 
using a Clean -box around the decrement, followed by a Clean over the entire field, but never 
cleaning deeper than about 2a. 

CLEANed maps of the SZ effect are useful for showing whether we have detected any re- 
solved structure in the decrement, and for overlaying on X-ray or optical images to check for 
positional coincidence. However, it is necessary to do all analysis for calculation of Hq in the aper- 
ture plane ( prainge et al. 200"!] ), since the noises are independent between visibilities, as opposed 
to the map plane in which the pixel noises are strongly correlated. 



3 RT OBSERVATIONS OF CL 0016+16 

Two array configurations were used: the most compact array, Cb, which provides short baselines 
at hour angles around zero, but in which the antennas shadow each other at extreme hour angles, 
and a more extended array, Cd, to provide unshadowed north-south baselines. Figure [I] shows the 
aperture-plane coverage. The asymmetry is due to the ~ 3° offset of the baselines from east-west. 
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The cluster was observed for the total of 46 days in the most compact array, and 32 days in the 
more extended array, between 1993 July 20 and 1995 January 10 with simultaneous observations 
of 0007+171 as a phase calibrator. The pointing centre was 00'^15™58^5 + 16°09'42" (B1950.0), 
which was also the phase centre. Previous long integrations with the RT have shown that it does 
not suffer from systematic offsets at the phase centre. 

A map of the long-baseline data detected the presence of two sources in the field. Both of these 
correspond with sources in the survey of Moffet & Birkinshaw( |1989| ) (hereafter MB). The only 
other source within 6' (the RT primary beam) of our pointing centre that was detected by MB is 
source MB 14, the extended cluster halo emission. This will be discussed further in section [431 . 
Fluxfitter was run over the two array configurations separately (because of the possibility of 
source variability as discussed in Grainge et al. ( |1996| )), and the results are shown in table 0. Once 
these sources had been subtracted we mapped and CLEANed the entire data set at two different 
resolutions (Figures and ^. These maps show clear detections of the SZ decrement and the high- 
resolution image has a resolved east-west extension. The position of the peak and the extent of the 
SZ decrement agree well with the X-ray emission. 



4 DETERMINING THE HUBBLE CONSTANT 

The method we use to measure Hq from RT cluster data has been described in detail by Grainge 



et al ( |2001D . As in that paper we make our models assuming an Einstein-de Sitter cosmology, and 



then calculate corrections for other cosmologies. Here we give details relevant to CI 0016-1-16. 



4.1 X-ray image 

We use the public Rosat PSPC image of CI 0016-1-16 (PI Hughes), with live time of 43157 s, 
to fit for the shape of the gas density distribution in the cluster and to determine a value of 
h^/Hq. We use the 0.5-2 keV band to minimise the effect of galactic emission and absorption; 
for the later we assume a column density of 5.6 x lO^'^cm"^. We fit to a two-dimensional p- 
model using Poisson statistics as the number of counts per pixel is typically small, and excluding 
regions containing X-ray point sources. We estimate the mean background level in a 15" pixel 
to be 4.04 X 10^^ counts s^^. We find best-fit parameters of P = 0.714, core radii of 40.8" and 
32.9" at a position angle of 53° and a central electron density of uq = 1.373 x lO^h^/^ mT^, 
where Hq = 100/zkm^^Mpc^^ assuming that the line-of-sight depth through the cluster is the ge- 
ometric mean of the two elliptical axes in the plane of the sky and an X-ray emission constant 
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of 3.80 ±0.27 X lO"*^^ counts s^^ from 1 of gas of electron density 1 m ^ at a luminosity 
distance of 1 Mpc. These are in good agreement with the two-dimensional models of Hughes & 
Birkinshaw ( |1998| ), who use a very similar fitting method, and also with those of Neumann & 
Bohringer ( |1997| ), who fit with a statistic after a Gaussian smoothing. 

Figure ^ shows the ROSAT PSPC image, our X-ray model, and the residuals from our fit. The 
reduced value of for this fit is 0.945. However, since the data are Poisson distributed with 
a relatively low mean for many pixels, approximating the distribution to be Gaussian is a poor 
assumption and therefore is not the correct statistic to use. We attempt to evaluate the goodness 
of fit by using the mean of the log likelihoods of each pixel, given the model, over the image. 
From our model we calculate 100 realisations of the image using the Poisson distribution with 
the appropriate mean at each pixel, and find the distribution of L to be L = —2.339 ±0.034. The 
corresponding value for the RoSAT image is L = —2.32. We therefore conclude that the fit is a 
good one. 

4.2 X-ray temperature 

We adopt the X-ray temperature calculated by Hughes & Birkenshaw ( |1998D of 7.55]']q gg keV. 



We assume the cluster to be isothermal at this temperature. The effect of this probably inaccurate 



assumption on our estimate of Hq is discussed in Maggi et al. ( [200 1[ ) 



4.3 Cluster halo emission 

Giovannini & Feretti ([2000[) have mapped the cluster halo source MB 14 with the VLA in B, C 



and D array at L-band (1.4 GHz). These measurements give access to short enough baselines to 
ensure that the full flux density has been found and has not been resolved out. They find that the 
total integrated flux from this extended source is 5.5 mJy. Assuming that the halo has a similar 
spectrum to that found in A665 dJones and Saunders 1996[ ) we estimate that the halo's flux density 



will have fallen to 55 ^Jy by 15 GHz. We model the halo emission as a Gaussian with FWHM 2' 
in order to estimate the extent to which we will resolve out the extended emission. We find that 
only 51 ± 9% of the total zero-spacing flux is detected by the RT shortest baseline of 870 X. 

The effect that halo emission will have on our calculation of Hq is to reduce the flux density 
of the detected SZ effect, which will lead to an overestimate of Hq by {Shaio/SszY ^- Estimating 
that the SZ decrement detected by the RT has been reduced by 30 ± 30 /jJy we apply a correction 
of —7 ± 7% to our value of Hq. 
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8 Keith Grainge et al. 
4.4 Ho from CI 0016+16 

Combining our X-ray model and SZ data in the aperture plane we find a maximum likelihood fit 
of//o = 69lj^kms-i Mpc"i,no= 1-141 x 10^ m^^ and a central decrement of 1077 ± 110 //K. 
Figure ^ shows the likelihood of the SZ data as a function of Hq. We find that the data fit very 
well the model predictions with a value of reduced = 1.047 and that the SZ effect is detected 
with good significance on many baselines (Figure |^). Including the dominant sources of error 
(X-ray temperature, X-ray fitting, SZ, ellipticity) ( [Grainge et al. 20()T| ) we find final estimates of 



Ho = km s^^ Mpc^^ for a = 1.0 cosmology; Hq = SAtfg km s^^ Mpc^^ for Qm = 0.3 

-23 
-18 



and = 0.7; and Hq = 76+^0 km s^^ Mpc^^ for = 0.2 and = 0.0 



5 DISCUSSION 



Reese at el ( pOOOD and Hughes and Birkinshaw ( [T998[ ) both use CI 0016+16 to calculate Hq. 



_i_ 1 9 _i_9 111 

They find 63_g km s^ Mpc^ (in combination with data from MS 0451.6-0305; assuming 
Q.M = 0.3, D.A = 0.7) and 47+j2 km s"^ Mpc"^ (assuming Q.m = 0.2, Qa = 0.0) respectively. 
Many of the errors that go into these estimates are identical to the ones that affect our determination 
(e.g. the error on the cluster temperature). Therefore the best way to check the consistency between 
these results is to examine the various determinations of the SZ central decrement. These are 
1242 ± 105 /jK (Reese et al.) and 1207 ± 190 fjK (Hughes and Birkinshaw), and so are somewhat 
higher than our value of 1077 ±110 /uK. Our result and that of Reese at el. differ by 1 . lo and so are 
not significantly discrepant, but we briefly discuss possible systematic differences between the two 
data sets. Firstly it is possible that the halo emission (section p3| ) has a less steep spectrum than the 
value we have adopted. In order to account fully for the difference we calculate that the spectrum 
of the extended emission between 1.4 and 15 GHz would have to be approximately =1.2 
(5 oc v^"). This is flatter than the flattest- spectrum halo emission seen in clusters between 1.4 and 
5 GHz of ttj 4 = 1 .3 ( Panisch 1982D , and we would expect the spectrum to steepen further between 



5 and 15 GHz due to spectral aging (dominated by inverse Compton losses). Furthermore such a 
source would also affect the measurement of Hughes and Birkinshaw at 18 GHz to almost the 
same extent as our observations at 15 GHz, but their value is in very good agreement with Reese et 
al. A second possibility is contamination from source MIO. Reese et al. do not remove this source 
from their data but do estimate the effect it could have by extrapolating from the 1 .4 GHz flux 
and find that its contribution is negligible. We note that there is some evidence that this source is 
variable with the result that estimating fluxes from observations at different epochs is insecure, and 
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that its spectral index between 1.4 GHz and 15 GHz is approximately 0.1, much flatter than the 
assumed value of 0.7. We estimate that the effect of MB 10 might be as much as 75 /jK at 30 GHz 
and that since this source lies in the negative sidelobes of the point spread function with respect to 
the cluster centre, its removal would reduce their calculated central decrement. 



6 CONCLUSIONS 



We have estimated the Hubble Constant by combining SZ and X-ray data from the cluster CI 0016+16 
and find Hq = 69+jg km s"^ Mpc"^ for a Qm = 1-0 cosmology; Hq = 84+ km s"^ Mpc"^ for 
Q.M = 0.3 and Q./^ = 0.7. This is in good agreement with other recent determinations from a va- 
riety of methods: the Hubble Key Project result of //q = 72 ± 8 ( preedman et al. 2000| ); from a 



low-redshift sample of SZ clusters Hq = 64+ |j ( [Mason et al. 200 1| ); from a high-redhift sample of 



SZ cluster s Hp = 5SZi dJones et al. 2001[); from gravitational lensing Hp = 59Zi =1= 15 (F assnacht 



et al. 1999) and Hq — 69+}q QBiggs et al. 1999D ; and from combining cluster velocities, cosmic 



microwave background and supernovae results Hq = T^^\q ( [Bridle et al. 200 1[ ). Although optical 
data of CI 0016-1-16 shows evidence that it has undergone a recent merger event ( |Smail et al. 1995| ) 
we are able to fit a good X-ray model to this cluster. This suggests that it should be possible to 
estimate the value of the acceleration parameter, qo, of the universe by combining SZ and X-ray 
observations of a complete sample of high redshift (z > 0.5) clusters. This would be a useful check 



of the result from observations of high redshift supernovae ( Perlmutter et al. 1999[ ). 
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Configuration Ae 1 Ae 2 Ae 3 Ae 4 Baselines Available 



Ca 
Cb 
Cc 
Cd 



12 
12 
16 
32 



9 
10 
12 
24 



7 

8 
8 

16 



4 
4 
4 
10 



2,3,3,4,5,5,7,8,9,12 

2,2,4,4,4,6,8,8,10,12 

4,4,4,4,8,8,8,12,12,16 



6,8,8,10,14,16,16,22,24,32 



Table 1. Compact configurations of tiie Ryle Telescope. The RT antennas are numbered 1-8 from east to west, with antennas 1-4 
being mounted on a railway track. The convention used for specifying the parking positions is to label them relative to the closest fixed 
antenna 5, in units of 9 m (= 450X). There are stations every 36 m along the track at positions 4, 8, 12 ... 128, and additional ones at 
7, 9 and 10. Antennas 6, 7 and 8 are at positions —128, —256 and —384, but were not used in these observations. 



© 0000 RAS, MNRAS 000, 000-000 



12 Keith Grainge et al. 

Aperture Plane 




© 0000 RAS, MNRAS 000, 000-000 



The SZ effect in CI 0016+16 

Table 2. Sources subtracted from CI 0016+16 visibility data. Positions are B1950. The errors in the fluxes are 39^jJy and 41^Jy in arrays Cb 
Cd respectively. 

RA Dec flux density flux density MB 

Cb array (ply) Cd array (jjjy) designation 
00 15 56.2 16 04 04 417 187 IVIB 15 

0015 49.3 16 1011 234 139 IVlB 10 
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Figure 2. High resolution RT map of tlie SZ effect in CI 0016+16 (contours) superposed onto the ROSAT PSPC image. The contours are -250, 
-200, -150, -100, -50 and 50 ^Jy. The restoring beam is 143 x 51" at an angle of 4.3°. The rms pixel noise on the map is 40 /jjy. The entire 
source-subtracted data set has been used to produce this map, with a Gaussian taper with 1 /e points at 2 Hk. 
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Figure 3. Low resolution RT map of the SZ effect in CI 0016+16 (contours) superposed onto the ROSAT PSPC image. The contours are -800, -720, 
-640, -560, -480, -400, -320, -240, -160, -80 and 80 /jjy. The restormg beam is 333 X 103" at an angle of 1.6°. The rms pixel noise on the map is 
90 /dy. Only data from the shortest RT basehnes (< 870 X) have been used in this map. 
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Figure 5. Binned real part of the source-subtracted visibilities (filled points with error bars) for CI 0016+16 observation against best fit model 
(open points joined by Une). 
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Figure 6. Likelihood plot for values of ffo bom the fit of SZ decrement to the X-ray derived cluster model. The width of the plot indicates only the 
errors from the visibility data. 
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